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In this work we show that the presence of a vector field on cosmological scales could explain 
the present phase of accelerated expansion of the universe. The proposed theory contains no di- 
mensional parameters nor potential terms and does not require unnatural initial conditions in the 
early universe, thus avoiding the so called cosmic coincidence problem. In addition, it fits the data 
from high-redshift supernovae with excellent precision, making definite predictions for cosmological 
parameters. Upcoming observations will be able to clearly discriminate this model from standard 
cosmology with cosmological constant. 
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Recent cosmological observations 0, 0, H, 0] indicate 
that the universe is undergoing a phase of accelerated 
expansion. The fact that the present rate of expansion 
is accelerating rather than decelerating poses one of the 
most important problems of modern cosmology. Indeed, 
in Standard Cosmology our universe starts expanding af- 
ter the Big Bang, but the attractive nature of gravity for 
ordinary matter and radiation necessarily slows down the 
expansion rate. In order to have acceleration, Einstein's 
equations require the universe to be dominated by some 
sort of non-ordinary energy (usually called dark energy) 
with the particular property of having negative pressure. 

Although its nature is unknown, a simple phenomeno- 
logical description in which dark energy is understood as 
a cosmological constant, i.e. a perfect fluid with equa- 
tion of state, pa — — Pa, where pa and p\ are the 
energy density and pressure respectively, seems to fit 
observations with very good precision (ACDM model). 
Thus, ACDM suggests that around a 70% of the energy 
density of the universe today would be in the form of 
dark energy, whereas the remaining 30% would be non- 
relativistic matter (the contribution from radiation and 
curvature being negligible) 0. 

However the fact that today matter and dark en- 
ergy have comparable contributions to the energy den- 



sity (both around (10 eV) in h 



1 units) turns 



out to be difficult to understand if dark energy is a true 
cosmological constant. Indeed, the energy density of a 
cosmological constant remains constant throughout the 
history of the universe, whereas those of the rest of com- 
ponents (matter or radiation) grow as we go back in time. 
Then the question arises as to whether it is a coincidence 
(or not) that they have comparable values today when 
they have differed by many orders of magnitude in the 
past. In addition, the cosmological constant exhibits an- 
other related problem. Its scale (around 10 -3 eV) is more 
than 30 orders of magnitude smaller than the scale of the 
other dimensional constant appearing in the gravitational 
equations, G = Mp 2 with Mp ~ 10 19 GeV, and it is also 
difficult to explain from particle or other known physics. 
In order to avoid these problems several models have 



been proposed in which dark energy is a dynamical com- 
ponent rather than a cosmological constant. Such models 
are usually based on cosmological scalar fields or modi- 
fications of Einstein's gravity 0, S, S S|- However, in 
order to have acceleration at the right time, they typ- 
ically introduce unnatural dimensional scales, resulting 
once again in fine tuning or coincidence problems. 

In this paper we consider a completely different type 
of dark energy model which is not based on scalar fields, 
but in the dynamics of a vector field. Unlike previous 
works [^, [n| , it is shown that vector fields can give rise 
to periods of acceleration even in the absence of poten- 
tial terms (see ll( for a general analysis). The existence 
of such solutions does not rely on the introduction of 
complicated functions of the fields and its derivatives, 
but can be obtained with the simplest kinetic terms, in- 
cluding two fields and two derivatives, so that the model 
does not contain dimensional parameters. Furthermore, 
we show that the required initial conditions for the vector 
fields are natural. 

Let us start by writing the action of our vector-tensor 
theory of gravity containing only two fields and two 
derivatives and without potential terms [l3 |: 



S = 




R» V A»A V 



(1) 



Notice that the theory contains no free parameters, the 
only dimensional scale being Newton's constant. The nu- 
merical factor in front of the vector kinetic terms can 
be fixed by the field normalization. Also notice that 
RfivA^A" can be written as a combination of derivative 
terms as V V A U - V^V^. 

The classical equations of motion derived from the ac- 
tion in |T]) are the Einstein's and vector field equations: 



Rfiu — -^Rg^v — 8TrG(Tfj, u 
nA M +R ia ,A u = 0, 



T ) 



(2) 
(3) 
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where T uv is the conserved energy-momentum tensor for 
matter and radiation and T^ v is the energy-momentum 
tensor coming from the vector field. In this work we 
shall solve these equations for the simplest isotropic and 
homogeneous flat cosmologies. Thus, we assume that 
the spatial components of the vector field vanish, so that 
A n = (A (t), 0, 0, 0) and, therefore, the space-time geom- 
etry will be given by the flat Robertson- Walker metric: 



ds 2 



dt 2 



a 2 {t)8ijdx 1 'dx\ 



For this metric reads: 



(4) 



the Hubble parameter Hq and the constant A$ during 
radiation (which fixes the total matter density Hjvf) are 
specified, the model is completely determined. In other 
words, this model contains the same number of parame- 
ters as ACDM, i.e. the minimum number of parameters 
of a cosmological model with dark energy. As seen from 
Fig. 1 the evolution of the universe ends at a finite time 
tend with a singularity in which with a end fi- 

nite, p DE — ► oo, pde — * -oo and A (t end ) = Mp/(4s/tt). 
This corresponds to a Type III singularity according to 
the classification in 11311. 



A + 3HA - 3 2H 2 + H A o = 



(5) 



where H — a/a is the Hubble parameter. 

Assuming that the universe has gone through radiation 
and matter phases in which the contribution from dark 
energy was negligible, we can easily solve these equations 
in those periods just taking H = p/t, with p = 1/2 for 
radiation and p = 2/3 for matter eras respectively, which 
is equivalent to assuming that a oc t p . In that case, the 
above equation has a growing and a decaying solution: 



A (t) = A+t a ++A-t° 



(6) 



with Aq constants of integration and a± = — (l±l)/4 in 
the radiation era, and a± = (— 3± v / 33)/6 in the matter 
era. 

On the other hand, the (00) component of Einstein's 
equations reads: 



H 2 



8ttG 



PA 



(7) 



with a = M, R and: 



Pa = 3 -H 2 A 2 +ZHA A n - l -A 2 . (8) 



Using the growing mode solution in ([5]), we obtain: 



PA = PAoa , 



(9) 



with n — —4 in the radiation era and n = (v 33 — 9)/2 ~ 
— 1.63 in the matter era. Thus the energy density of the 
vector field starts scaling as radiation at early times, so 
that pa/ Pr — const. However when the universe enters 
its matter era, pa starts growing relative to pu eventu- 
ally overcoming it at some point, in which the dark en- 
ergy vector field would become the dominant component. 
From that point on, we cannot obtain analytic solutions 
to the field equations. In Fig. 1 we show the numerical 
solution to the exact equations, which confirms our ana- 
lytical estimates in the radiation and matter eras. Notice 
that since An is essentially constant during the radiation 
era, solutions do not depend on the precise initial time 
at which we specify it. Thus, once the present value of 



p/p 10 30 




FIG. 1: Evolution of energy densities for the best fit model. 
Dashed (red) for radiation, dotted (green) for matter and solid 
(blue) for vector dark energy. We show also for comparison 
the cosmological constant density in dashed-dotted line. We 
see the rapid growth of the vector dark energy contribution 
at late times approaching the final singularity. 



We can also calculate the effective equation of state for 
dark energy as: 



wde 



PA 
PA 



-3 ( %H l 



h)a 2 + ha a -IA 2 



(10) 



\R^A 2 + 3HA An - \A 2 
Again, using the approximate solutions in ([S]), we obtain; 



wde 



3 

V33-15 



radiation era 
-0.457 matter era 



(11) 



After dark energy starts dominating, the equation of 
state abruptly falls towards woe — > — 00 as the universe 
approaches t end . As shown in Fig. 2 the equation of state 
can cross the so called phantom divide, so that we can 
have u>de{z = 0) < — 1. In Fig. 3, we show the evolution 
of the Aq component. 

In order to confront the predictions of the model with 
observations of high-redshift type la supernovae, we have 
calculated the distance modulus as a function of redshift. 
Comparing pth(z) with its observational value in a given 
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FIG. 2: Evolution of dark energy equation of state for the best 
fit model. The lower panel shows the la confidence interval. 
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TABLE I: Best fit parameters with lcr intervals for the vector 
model (VCDM) and the cosmological constant model (ACDM) for 
the Gold (157 SNe) and SNLS (115 SNe) data sets. w denotes the 
present equation of state of dark energy, ylo is the constant value of 
the vector field component during radiation. 2y is the deceleration- 
aceleration transition redshift. to is the age of the universe in units 
of the present Hubble time. t enl i is the duration of the universe in 
the same units. 



FIG. 3: Evolution of Ao in Mp units for the best fit model. 



data set will enable us to carry out a % 2 statistical anal- 
ysis. For this purpose, we have considered two sets of 
supernovae: the Gold set 14[, containing 157 points with 
2 < 1.7, and the more recent SNLS data set [Hj|, com- 
prising 115 supernovae but with lower redshifts (z < 1). 

In Table 1 we show the results for the best fit together 
with its corresponding lcr intervals for the two data sets. 
We also show for comparison the results for a standard 
ACDM model. We see that the vector model (VCDM) 
fits the data considerably better than ACDM (at more 
than 2a) in the Gold set, whereas the situation is reversed 
in the SNLS set. This is just a reflection of the well- 
known 2a tension [l6| between the two data sets. The 
best fit parameters for the VCDM model are identical for 
the two data sets with small differences in the confidence 
intervals. Compared with ACDM, we see that VCDM fa- 
vors a younger universe (in Hq 1 units) with larger matter 
density. In addition, the deceleration-aceleration transi- 
tion takes place at a lower redshift in the VCDM case. 
Another important difference arises in the present value 
of the equation of state with wq = — 3.53^Q gy which 
clearly excludes the cosmological constant value —1. Fu- 



ture surveys [17[ are expected to be able to measure Wq 
at the few percent level and therefore could discriminate 
between the two models. 

We have also compared with other parametrizations for 
the dark energy equation of state. Thus for instance, tak- 
ing wde(z) — wo + wiz(l + z)" 1 [18}, we find \ 2 = 173.5 
for the Gold set. Since this is a three-parameter fit, in 
order to compare with the one-parameter fits of VCDM 
or ACDM, we use the reduced chi-squared: x^/d.o.f — 
1.108 for VCDM, x 2 /d.o.f = 1.127 for the (w ,wi) 
parametrization and x^/d.o.f = 1.135 for ACDM. As 
a matter of fact, to our knowledge best, VCDM provides 
the best fit to date for the Gold data set, since the oscil- 
latory four-parameter model previously reported in [l9[ 
still has x 2 /d.o.f = 1.115. 

The evolution of dark energy for the best-fit model 
is plotted in Figs. 1-3. We see that unlike the cos- 
mological constant case, throughout the radiation era 
Pde/pb. ~ 10~ 6 . Notice that although the onset of cos- 
mic acceleration depends on the value of Aq during that 
era, for the best-fit Ao = 3.71 x 10 -4 Mp, which is rela- 
tively close to the Planck scale and could arise naturally 
in the early universe without the need of introducing ex- 
tremely small parameters. 
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When comparing the parameters obtained from SN 
la (Table 1) with predictions coming from CMB 
anisotropics or baryon acoustic oscillations [2(J, it is im- 
portant to keep in mind that such predictions are ob- 
tained after a data process which involves the use of a 
particular model for dark energy, which in most cases is 
ACDM, as a fiducial model. This is a good approxima- 
tion for models with (nearly) constant equation of state 
[20l | . but could not be a priori justified in our case since 
wde(z) has a strong redshift dependence [21 1. 

In this work we have only considered the time compo- 
nent of the vector field. The presence of spatial compo- 
nents could, in principle, have adverse effects. However, 
we have found that the energy density of the spatial part 
decays as a -8 during radiation and matter eras, i.e. much 
faster than the temporal contribution, so that it will not 
dominate at late times. On the other hand, we have 
calculated the evolution of (pn — p±)/p, where pu is the 
pressure along the direction of the spatial component and 
p± is the transverse pressure, and we have found that this 
quantity decays very fast during the matter and radia- 
tion eras. Accordingly, we do not expect the generation 
of large anisotropics. 

So far we have only considered the homogeneous 
model. In order to study the model stability we have 
considered the evolution of metric and vector field per- 
turbations. Thus, we obtain the dispersion relation and 
the propagation speed of scalar, vector and tensor modes. 
For all of them we obtain v = (1 — 167tGAq) -1 / 2 which 
is real throughout the universe evolution, since the value 
Aq = (WttG)^ 1 exactly corresponds to that at the final 
singularity. Therefore the model does not exhibit expo- 
nential instabilities. As shown in [2^], the fact that the 
propagation speed is faster than c does not necessarily 
imply inconsistencies with causality. We have also con- 
sidered the evolution of scalar perturbations in the vector 
field generated by scalar metric perturbations during the 
matter and radiation eras, and found that the energy den- 
sity contrast dp a I PA is constant on super-Hubble scales, 
whereas it oscillates with growing amplitude as a 2 in the 
radiation era and as ~ a 3 in the matter era for sub- 
Hubble scales. Therefore again, we do not find exponen- 
tially growing modes. 

The model proposed in this work can be considered 
as an effective description of dark energy on cosmologi- 
cal scales. Extending the applicability range to smaller 
scales requires consistency with local gravity tests. In- 
deed, we can see that for the model in JT]), the static 
post-Newtonian parameters agree with those of General 



Relativity (12|, i.e. 7 = j3 = 1. For the parameters as- 
sociated to preferred frame effects we get: a\ = and 
a-2 = 8itAq/Mp where Aq is the norm of the vector 
field at the solar system scale. Current limits ct2 ~ 10 -4 
(or a2 < 1CP 7 for static vector fields during solar system 
formation) then impose a bound Aq< 10~ 5 (10 -8 ) Mf>, 
which could conflict with the model predictions, since 



the present (Solar system formation) values on cosmo- 
logical scales are: 1.3 x 1CT 1 M P (7.5 x 1CT 2 M P ). How- 
ever, notice that the cosmological values do not need to 
agree with those at lower scales. The latter will be deter- 
mined by the mechanism that generated this field in the 
early universe characterized by its primordial spectrum 
of perturbations, and the subsequent evolution in the for- 
mation of the galaxy and Solar system. Concerning the 
potential presence of quantum instabilities in the model, 
in (23l | the condition in order to ensure positive norm 
Hilbert space are obtained. With our Riemann tensor 
sign convention, we see that such a condition is indeed 
satisfied in the model ([1]). 

In conclusion, the results of this work show that vector 
theories offer an accurate phenomenological description 
of dark energy in which fine tuning problems could be 
avoided. 
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